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Abstract. We consider a question of finding a periodic solution for the planar 
Newtonian N-body problem with equal masses, where each body is travelling 
along the same closed path. We provide a computer assisted proof for the 
following facts: local uniqueness and convexity of Chenciner and Montgomery 
Eight, an existence (and local uniqueness) for Gerver's SuperEight for 4-bodies 
and a doubly symmetric linear chain for 6-bodies. 



1. Introduction 



In this paper we consider the problem of finding periodic solution to the iV-body 
problem in which all N masses travel along a fixed curve in the plane. The A^-body 
problem with N equal unit masses is given by a differential equation 



(1.1) 



, N, n 



= ||(7i — (7j||. The gravitational constant is taken 



where qi £ M", i — \ 
equal to 1 . 

We consider planar case (n=2) only, we set qi = {xi, yi), qi — {vi, Ui). Using this 
we can express (jl.lf) by: 



(1.2) 



Ui 



Recently, this problem received a lot attention in literature see |M] ICMI ICGMSI 
EHISSIHIIMH and papers cited there. 

By a simple choreography |S1[IS2| we mean a collision- free solution of the N-body 
problem in which all masses move on the same curve with a constant phase shift. 
This means that there exists q : K — > a T-periodic function of time, such 
that the position of fc-th body {k — 0, . . . , — 1) is given by qkit) = q{t + k-^) 
and (goj 911 • ■ • j qN-i) is solution of the A'^-body problem. The simplest choreogra- 
phies are Langrange solutions in which the bodies are located at the vertices of a 
regular N-gon and move with constant angular velocity. Another simple choreog- 
raphy, a figure eight curve (see Figure^, was found numerically by C.Moore fM|. 
A. Chenciner and R.Montgomery jCMj gave a rigorous existence proof of the Eight 
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in 2000. In December 1999, J.Gerver found orbit for iV = 4 called 'Super-Eight' 
(Figure EJ. After that C.Simo found a lot of simple choreographies with different 
shape and number of bodies ranging form 4 to several hundreds (see [SJ ISII ICGMS| 
for pictures, animations and more details). 

Up to now the only choreographies whose existence has been rigorously estab- 
lished MR are Lagrange solutions and the Eight solution. While the Lagrange 
solution is given analytically, the existence of Eight was proven in (CM) using vari- 
ational arguments and there is still a lot of open questions about it jChI IMR| . For 
example uniqueness (up to obvious symmetries and rescaling ) and convexity of the 
lobes in Eight. In SectionOwe give a computer assisted proof of existence of Eight, 
its local uniqueness and convexity of the lobes. 

In Sect ions 01 island El we concentrate on choreographies called doubly symmetric 
linear chains. They are symmetric with respect to the x and y coordinate axes 
and all points of the curve self-intersection are on the x axis. We give a computer 
assisted proof of an existence (and local uniqueness) of doubly symmetric linear 
chains for four (Gerver SuperEight) and six bodies. 

Proofs given in this paper are computer assisted. By this we mean that we 
use computer to provide rigorous bounds for solutions of The problem of 

proving an existence of a choreography is reduced to finding a zero for a suitable 
function. For this purpose we use an interval Newton method and a Krawczyk 
method |XllKllMollN| . (see Sectionj^) which are apparently rather unknown outside 
the interval arithmetic community. To integrate equations we use a C^-Lohner 
algorithm |ZLo| . All computations were performed on AMD Athlon 1700XP with 
256 MB DDRA M memory, with Windows 98SE operating system. We used CAPD 
package |Capd| and Borland C-I--I- 5.02 compiler. A total computation time for 6- 
bodies was under 90 seconds and considerably smaller for the Eight and SuperEight 
solutions - see Section [3 for more details. 

2. Two ZERO FINDING METHODS 

The main technical tool used in this paper in order to establish an existence of 
solutions of equations of the form f{x) = is on interval Newton m,ethod \K\ IMol IN] 
and a Krawczyk method^^ El IH] ■ The interval Newton method was used to prove 
an existence of the Eight (Figure and Gerver orbit (Figure El . The Krawczyk 
method for Gerver orbit and an orbit with 6 bodies in a linear chain (Figure ISJ. 

2.1. Notation. In the application of interval arithmetics to rigorous verification 
of theorems single valued objects, like numbers, vectors, matrices etc are in the 
formulas replaced by sets containing sure bounds for them. In the sequel, we will 
not use any special notation for single valued object and sets. For a set S by [S] 
we denote the interval hull of S, i.e. the smallest product of intervals containing S. 
For a set which is an interval set (i.e. can be represented as a product of intervals) 
we will also use square brackets to stress its interval nature. For any interval set 
[S] by mid{[S]) we denote a center point of [S]. For any interval [a,b] we define a 
diameter by diam[a, b] := b — a. For an interval vector (matrix) S = [S] by diamS 
we denote a vector (matrix) of diameters of each components. 

2.2. Interval Newton method. 

Theorem 2.1. Let F : M" M" be a function. Let [X] = nf^Ja,,6,], 

Ui < bi. Assume the interval hull of DF([X]), denoted here by [DF{[X])], is 
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invertible. Let x £ X and we define 

(2.1) N{x,[X])^-[DFi[X])]-'Fix) + x 
Then 

0. if Xi,X2 G [X] and F(xi) — F{x2), then xi — xi 

1. if N{x, [X]) C [X], then 3\x* G [X] such that F{x*) = 

2. if xi e [X] and F{xi) = 0, then Xi G N{x, [X]) 

3. if N{x, [X]) n [X] = 0, then F{x) ^ /or all x G [X] 

2.3. Krawczyk method. We assume that: 

• F : R" — > R" is a function, 

• [X] C R" is an interval set, 

• xe[x] 

• C G R"^" is a linear isomorphism. 
The Krawczyk El IH] operator is given by 

(2.2) Kix, [X],F) -.^x- CF{x) + [Id - C[DF{[X])]){[X] - x). 

Theorem 2.2. 1. If x* G [X] and F{x*) = 0, then x* G [X],F). 

2. If K(x,[X],F) C mt[X], then there exists in [X] exactly one solution of 
equation F(x) = 0. 

3. IfK{x, [XlF) n [X] = 0, then F{x) ^ /or all x G [X] 

2.4. Algorithm for Newton and Krawczyk method. Theorems 12.11 and 12.21 

can be used as a basis for an algorithm for rigorous enclosing for solution of equation 
F{x) = 0. Let T{x, [X]) = N{x, [X]) if we are using the interval Newton method 
and T{x, [X]) = K{x, [X], F) in case of Krawczyk method. 

First, we need to have a good guess for x* G R". For this purpose we use 
a nonrigorous Newton method to obtain x = {xi,X2, ■ ■ ■ ,Xn)- Then we choose 
interval set [X] which contains x and perform the following algorithm: 
Step 1. Compute T{x, [X]). 
Step 2. If T(a;, [X]) C [X], then return success. 

Step 3. If X n T{x, [X]) = 0, then return fail. There are no zeroes of F in [X]. 
Step 4. If [X] C T(x, [X]), then modify computation parameters (for example: a 

time step, the order of Taylor method, size of [X]). Go to Step 1. 
Step 5. Define a new [X] by [X] := [X] n T{x, [X]) and a new x hy x := mid{[X]), 

then go to Step 1. 

In practical computation it is convenient to define a maximum number of itera- 
tion allowed and return fail if the actual iteration count is larger. 

Observe that the third assertion in both Theorems 12.11 and 12.21 can be used to 
exclude an existence of zero of F. This has been used by Galias |G1[ IG2| to find 
all periodic orbits up to a given period for Henon map and Ikeda map. It seems 
possible to obtain similar results for N-body problem in the future. 

3. The Eight - existence, local uniqueness and convexity 

Existence of the Eight has been shown in jCMj using mixture of symmetry and 
variational arguments. Here we give an another existence proof and in addition 
we obtain local uniqueness and the convexity of each lobe of the Eight. We follow 
|CM| in the use of the symmetry, but other component of the proof is different - 
we use an interval Newton method discussed in Sectional 
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In notation, but only for the Eight, we follow |CMj . 
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Figure 1. The Eight - the initial position 

Let T be any positive real number. We define action of the Klein group Z2 x Z2 
on R/TZ and on as follows: if a and t are generators, 
T T 

(3.1) a{t) + —, Tit) = + 2" ' v) = (~^' y)' '^(^' y) = ~y)- 

For loop q : (R/TZ) — > R^ and i — 1,2,3 we define the position for i-th body 

by 

(3.2) q^{t)^q(^+(Z-i).^^ 

Here qi is just position of i-th body. 

The following theorem without the uniqueness part was proved in |UM| . 

Theorem 3.1. There exists an "eight" -shaped planar loop q : (R/TZ) — > R^ with 
the following properties: 

(1) for each t, 

qi{t) + q2{t) + qsit) = 0; 

(2) q is invariant with respect to the action 0/Z2 x Z2 on R/TZ and on R^; 

q o a{t) = a o qit) and q o T(t) —to q{t); 

(3) the loop X : R/TZ — > R^ defined by 

x{t) = {qi{t),q2{t),q3{t)) 
is T-periodic solution of the planar three-body problem with equal masses. 











/ 1 " 
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Figure 2. The Eight - the final position 
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Moreover, q is locally unique (up to obvious symmetries and reseating). 

Remark 3.2. If conditions (1),(2),(3) are satisfied then 
a. qi{t) + q2{t) + q-?,{t) = for each t, 
h. qo a[t) — a o q(t) and q o T{t) = a o q(t) for each t, 

c. (73(0) = (0,0) and 93(0) = -2gi(0) 

d. gi(0) ^ -92(0) and gi(0) = 92(0), 

e. qi{T/12) is on the X axis and gi(T/12) is orthogonal to the X axis, 

f. qsiT/U) = T o (72(T/12) and 93(^/12) = a o g2(r/12), 

The following lemma describes the symmetry reduction for the Eight 

Lemma 3.3. Assume that q : [0,T] > is a solution of the three body problem, 

such that q = ('?i,92,93) satisfies conditions (c),(d),(e),(f) in Remark then 

exists q : {K./T'L) > satisfying condition (1),(2),(3) in Theorem with 

T = 12f. 

Proof: We define 

93 (i) forie[0,f] 
(3.3) q{t) = <( T o 92(2f - t) for t e [f, 2f] 



aoqi{t- 2T) for t e [2T, ST] 



and 



(3.4) q{t) = 



q{t) for t e [0, 3f ] 

Toq{T-\t)) forte[3f,6f] 

aoq{a-^t)) fortG[6f,9f] 

aoToq{T-^ oa~^{t)) for t £ [9f , 12f ] 

Let /i and /2 be two solutions of II. II on intervals [ii, i2] and [^2,^3] respectively. 
If fi{t2) = 72(^2) and /f (^2) = /2^(^2) then / = {/i,/2} is solution on interval 
[ti, t2]- To show that q{t) is a solution it is enough to show that "pieces fit together 
smoothly" . 

For t^f from Of we have 

(3.5) UT) = <73(J)=T092(J)=T092(f)=T092(2f-f) 

(3.6) 9V(f) = ao4+(f) 
For t = 2f from Eld wc have 

(3.7) 9(2f) = r 092(0) = a 091(0) 

(3.8) tim = ao^-{f)=ao^+{2f)=^+{2f) 

From (|3.5I3.8|I it follows that 9(t) is smooth curve for t € [0, 3T]. We construct 
q(t) from that curve. 
For t — 3T from l3.2l e we obtain 

9(3T) = (7 o 9i(T) = (T o T o 9i(T') = cr o T o 9(r^"'"(3T)) 

9-(3f) - ao9\(f)-9\(r)-9+(3f) 

Using c for t — 6T we infer that 

9(6f) - ro9(r-i(6f)) 093(0) = (7093(0) = ao9>-i(6f)) 

9-(6f) = ao9\(f) = 9+(6r) 
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Now we show that q{t) is a closed curve. 

q{12f) ^aoTo q{T-^ o a^'^iuf)) = cr o r o q{0) = q{0) = q{0) 

r (12f ) ^-[aoTo q{T-^ o a-\\2f))] ^ (j(0) = <7+(0) 

Hence it is easy to see that q{t) can be extended to a T-periodic curve, such that 
x{t) = {qi{t),q2{t),q3{t)) is a T-periodic solution of the three body problem. 

Condition (2) in Theorem |0] follows easily from definition of q{t) and properties 
of cr and r (for T-periodic orbit ^-^(t) a{t) and T-i(i) = r(i)). For t g [0,3f] 
we have 

(3.9) qicr{t)) = a o q{a-^ o a{t)) ^ a o q(t) = a o q{t) 

(3.10) qiT{t)) = Tog(T"l or(t)) = rog(t) rog(<) 

For t e [Sf, 6f ] using ifTTUIl we obtain 

(3.11) q{(T{t)) = (7 o T o q{T^^ o a^^ o a{t')) = a o t o q{T^^{t)) = a o q(t) 

(3.12) q{T{t)) q{T{t)) =Toq{To T{t)) - r o g(i) 

We omit other two cases, because the proof is very similar. 

Because there aren't any external force so the center of mass can only move with 
a constant velocity. But from 13. 21 c and 13. 21 d we have gi(0) -I- 92(0) -I- 93(0) = and 
(71(0) -I- (72(0) + 93(0) ~ 0. Hence we obtain condition (1) in Theorem l3.1l □ 




Figure 3. Rotated Eight - initial position 




Figure 4. Rotated Eight - final position 
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Hence to prove Theorem 13.11 it is enough to show that there exists a locaUy 
unique (up to obvious degeneracies) function satisfying assumptions of Lemma 13.31 
For this end we rewrite these assumptions as a zero finding problem to which we 
apply an interval Newton method in the reduced space. 

Our original phase space is 12 dimensional, the state of bodies is given by 
(qi, (72, 93, 91, 92, (/a)- The center of mass is fixed at the origin. Hence one body's 
position and velocity is determined by other two bodies. We start from a coUinear 
position with the third body at the origin and with equal velocities of the first and 
the second body (seeE)). Hence it is enough to know the position and the velocity 
of the first body to reconstruct initial condition of other bodies. Moreover, if we 
have one solution we could get another solution by a suitable rotation (both have 
the same shape). To remove this degeneracy we place the first body on the X axis 
(we will make computation for rotated Eight, see Figure |^. In addition we fix 
the size of trajectory. This fixes the period of the solution. But if we know one 
periodic solution then by Kepler law we may obtain a solution of any period just 
by rescaling. Finally, we set 91(0) = (1, 0), hence initial conditions are fixed by the 
velocity of the first body. 

The reduced space for Eight is two dimensional and is parameterized by velocity 
of second body. We define map from the reduced space to full phase space E : 
M? — > M}"^ , which expands velocity of first body, given by (f,u), to the initial 
conditions (2:1, yi, 2:2, 1/2, 2:3, j/3, wi, ui, W2, "2, ^^3, U3) for equation (|l.ll) 

E{v,u) = (l,0,-l,0,0,0,t;,M,'i;,u, -2u,-2u). 

For each such initial condition exists a solution of the tree body problem defined 
on some interval. To each initial configuration, following that solution, we asso- 
ciate, if it exists, a configuration in which for the first body the position vector is 
orthogonal to its velocity vector for the first time. This defines the Poincare map 

p -.M}^ — > M}^. 

Now, we define map R : M^^ — > R^, by 

^(91,92,93,91,92,93) = (II92 - 9i|P - II93 - 9i|P, (92 - 93) X 91), 
where by x we denote vector product, and map $ : M? — > by 

^ = Ro P o E. 

Remark 3.4. If i?(gi, 52, 93, 9i, 92, 93) = 0, then \\q2 - qi\\ = II93 - 9i|| and if in 
addition yi = but xi ^ 0, then y2 ~ ys- The bodies are then located as in Fig. ^ 
and conditions (e) and (f) in Rem. 13.21 are satisfied in a suitably rotated coordinate 
frame. 

The following lemma, which is crucial for the proof of Theorem 13.11 is obtained 
with computer assistance. 

Lemma 3.5. There exists a locally unique {v,u) G that <^{v,u) — (0,0). 

Proof: We use an interval Newton method fTheorcm I2.1|l . First we come close 
to zero, starting with some rough initial condition (for example from using a 
non-rigorous Newton method. Once we have a good candidate xq — {vq, uq), we set 
[X] = [vo — S,vo + 6] X [uo — 6,uo + S] and compute rigorously ^{xq) and ^^g^^- For 
this purpose we use C^-Lohner algorithm described in (ZCo) . In this computation 
we used a time step /i = 0.01 and the order r — 7. 
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Result 



([-2.107029e-06, -2.106467e-06],[2.974991e-06, 2.976034e-06]) 



diam ^{xq) 



(5.625889e-10, 1.042962e-09) 



dx 



[17.622624,17.643043] [1.809772,1.827325] 
[-24.868548, -24.848432] [-10.056629, -10.039221] 



N{xo,[X]) 



([0.347116886243943,0.347116889993313], 
[0.532724941587373,0.532724949187495]) 



diam Njxn, [X]) 



(3.749369e-09,7.600121e-09) 



Table 1. The data from the proof of Lemma [3. 51 



It turns out that the assumption of assertion 1 in Theorem 12 . II holds for 
xo = (0.347116768716,0.532724944657) and S = lO"*'. Numerical data from this 
computation are listed in Tabled 

Moreover, from Theorem 12. II we know that this zero is unique in the set X. □ 
Proof of Theorem 13.11 From Lemma f3. 51 there exists an initial condition {v,u) 
in reduced space that $(u, u) = (0, 0). Hence there exists solution q{t) of the three 
body problem in some interval [0, T] that in < = bodies are in coUinear position 
with third body in the origin and in i = T they are in isosceles configuration 
(ri2 = fia). By rotation we could get solution q{t) that in t = T first body is on 
X axis. Obviously we still start in coUinear position and end in isosceles. Now 
from Remark l3.4l gft') satisfy all condition in Lemma 13.31 and hence there exist q{t) 
satisfying condition (1),(2),(3). 

Local uniqueness follows from local uniqueness in Lemma 13.51 □ 

3.1. Convexity of the Eight. 

Theorem 3.6. Each lobe of the Eight is convex. 

Proof: For the proof it is enough to show that the only inflection point on the 
curve q{t) is the origin. 

From Lemma 13.51 we get set [X] which includes initial condition for the Eight 
in reduced space. We expand [X] to the full space, we set X = E{[X]). In the 
coordinate frame in which the Eight looks like in Fig. ^ and symmetries a and r 
are reflections with respect to coordinate axes fsee l3.1|l we see immediately that the 

symmetry properties of the Eight imply that at the origin 1^ = and 1^ = 0. 
To prove convexity of the Eight we follow rigorously the trajectory of set X and 
show that the only point in which = and = is the origin. The same 
is true if we rotate the Eight, as in Figure |2| to the coordinate system in which we 
performed the actual computations and in which we will work for the remainder of 
this proof. 

Let qi(t) = {xi{t),yi{t)) be position of i-th body . If ^ ^[tk-i,tk] then, 
from the implicit function theorem, we can write yi as a function of Xi for interval 
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between Xi{tk-i) and Xi{tk). Otherwise we try to represent Xi as a function of yi. 
For small enough time steps at least one of this representations is always possible 
for the Eight (this is really verified during the rigorous computations) . Now we 
derive the formulas for the derivatives of yi{t{xi)) with respect to Xi (for derivatives 
of Xi{t{yi)) we just exchange x and y variables in formulas). 

dyi{t{xi)) _ % / dx. 



dxi dt \ dt 

d^yi{t{xi)) f d^yi d'^Xi dyi\ / dx. 



dxf \ dt^ dt^ dx, J \dt 

d-^yi{t{xi)) 0*3 dt ~dW dt '^^\~dt^J dx, ^^tT df^ \dtT dxf ) 

^ ST 

From this equations we need to know the derivatives Xi and yi with respect to time. 
We obtain them easily by differentiation of Hl.l|) with respect to time. In fact this 
is done during each step of C^-Lohner algorithm. 

By (p{t, xq) we denote the state of bodies (the position and the velocity) at time 
t with the initial condition xq a,t t = 0. Let hk be length of k-th time step, = 

/ii H \-hk- the time after k steps, [q''] = [q'[] x [q^] x [q^] x [qf[] x x [gf ] C M.^"^ 

be product of intervals, such that ip{tk,X) C [q^] and [Q'^] = [Q^] x [Q2] x [Q^] x 
[Qi] X [Q2] X [Q3] C where [Qf] = [X^] x [F^], be product of intervals, such 
that ip{[tk-i,tk],X) C [Q''] Note that [Q''] can be seen as an interval enclosure for 
whole trajectory between [g*^^^] and [q''], which is computed during each step of 
C^-Lohner algorithm. 

For fc-th step and i-th body (except first step for third body) we check if at least 
one of the following conditions is true 



(3.13) o^-^[Xf]andO^— ^^[X 



^[Xt] andO^^' 



(3.14) o^^[r^]andO^^[r^]. 

For first step for third body (starting in the origin) we check if one of following 
conditions is satisfied 

(3.15) i ^[Xi] and £ ^[X^] and ^ ^[X^] 

(3.16) i ^[Yi] and £ ^[Yi] and ^ ^[Yi]. 

To verify above conditions we follow rigorously using C^-Lohner algorithm the 
trajectory of set X until reaching the section described in proof of Lemma [3. 51 and 
for each time step and each body we check suitable condition. It turns out that 
for each time step and each body these conditions were satisfied. This finishes the 
proof. □ 

3.2. Some numerical data from the convexity of Eight computation. The 

parameters of the methods were: a time step h = 0.01, order r = 7. We needed 53 
steps to cross the section, below we show data for characteristic cases in order to 
show that it was really quite easy to verify, with algorithms we used. 
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^ixh 


1 
2 
3 


[0.334402,0.347118] 
[0.347116,0.360049] 
[-0.695034,-0.69385] 


[15.3592,17.9897] 
[-16.5013,-14.111] 
[-0.0682713,0.269952] 


[136.616,219.114] 
[119.951,192.562] 
[-30.969,-26.3718] 



Table 2. The data from the proof of Theorem 13.61 for step 1, 
neighborhood of deflection point. 



i 


dxi r v2l 




d'''xi r v2l 
^lX^\ 


1 

2 
3 


[0.32222,0.334722] 
[0.35972,0.372882] 
[-0.695669,-0.69444] 


[16.7085,19.6225] 
[-15.2203,-13.0177] 
[0.126359,0.472046] 


[155.007,250.021] 
[105.778,171.191] 
[-30.9533,-26.1203] 



Table 3. The data from the proof of Theorem 13.61 for step 2. 



i 


^[Xf] 




^ixf] 


1 
2 
3 


[-0.00287209,0.00468403] 
[0.904939,0.922079] 
[-0.919428,-0.909617] 


[-2.55715,-2.16831] 
[2.56371,3.03259] 


[4.35197,10.6201] 
[-3.51203,3.48564] 


i 


9yi [^-21 

at i 


d^Vi ry2i 


'a^y^i J 


1 


[0.480975,0.48288] 


[-3.24824,-3.11737] 


[-2.98453,-0.860616] 



Table 4. The data from the proof of Theorem 13. 61 for step 37. 



Step 1. We start in colhnear position with third body in the origin. So we are in an 
inflection point. Numerical data for this case are given in Table |21 We see that for 
third body we have e ^^[X/], but this derivative is monotonic (^^[X/] < 0), 

hence there can be only one zero of it in interval [^3]. But we know that one zero 
is at (0,0) G [Xl]. 

Steps 2-36 and 38-53 In this case all second derivatives do not contain 0. TableEl 
contains data obtained in second step. 

Step 37. In this case we cannot represent yi as function of xi, hence we interchange 
variables and represent xi as function of yi and check condition l3. 14l instead of 13. 131 
The first body is in the rightmost position. Table 01 contains the derivatives for this 
case. 

4. Doubly symmetric choreographies with even number of bodies 

4.1. Symmetries. Many of the choreographies found by Simo |S1| have at least 
one symmetry. Just as in case of the Eight this is not only a symmetry of trajectory 
image, but also of how the curve describing the trajectory is parameterized with 
time. 
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Figure 5. Linear chain with 8 bodies. 



In this section we introduce a notation for symmetries which wiU be used till the 
end of this paper. By Sx, Sy, So we want to denote a symmetry with respect to 
the X axis, the y axis and to the origin. These spatial symmetries act also on time 
variable parameterizing curves. 

Let T be any positive real number. We define actions of Sx, Sy and 5*0 on M/TZ 
and on as follows: 

Sx{t)^-t+^, Sx(x,y) ^ {-x,y), 
Sy{t) = -t, Sy{x, y) = (x, -y), 

So{t)=t+'^, So{x,y) = {-x,~y). 

It follows from this definition that Sq = Sx ° Sy — Sy o Sx- 

Let q{t) : R/TZ — > R^. We say that q{t) is invariant (equivariant) with respect 
to the action of S if S{q{t)) ~ q{S{t)) for all t. Moreover we assume that q{t) is 
function. If q{t) is invariant with respect to Sx (resp. 5*^) then q(Sx{t)) = Sy{q(t)) 
(resp. q{Sy{t)) — Sx{q{t))). Hence Sq invariance implies that q{So{t)) = So{q{t)). 

From now on we will enumerate bodies starting from 0. We set qi{t) = q (t + ^ ■ z) 
fori = 0,...,iV-l. 

4.2. Doubly symmetric choreographies with even number of bodies. We 

will consider only cases with even number of bodies. 

Let T = TV • f > 0. We search for a function q{t) : R/TZ — > R^ which has 
following properties: 

PI. for each t the origin is center of mass, 

JV-l 

(4.1) ^g(< + T.z)=0, 

1=0 

P2. q{t) is invariant with respect to 

(a) Sx i.e. qiSxit)) = Sxiqit)), 

(b) ^0 i.e. qiSoit)) = Soiqit)), 

P3. function x — (go(i), 92(^)1 ^w-i) where qi{t) = q(t + T-i^ for i = 
0, . . . , N — 1, is a T-periodic solution of the Newtonian N-body problem 

Observe that condition P3 says that bodies trace each other with a constant 
phase shift. 

The following Lemma, which is analogous to Lemma 13.31 gives necessary and 
sufficient conditions for an existence of a choreography satisfying P1,P2 and P3. 

Lemma 4.1. Let N = 2n be a number of bodies. There exists a function q(t) with 
properties PI, P2, P3 if and only if there are functions qi : [0,T/2] — > for 
i = 0, l,...,iV — 1 such as: 

(1) go(0) = (0,2/0) for some j/o ^ (qo{0) is on Y axis), 
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(2) At time to = for i = 0,1, ... , N/2 - 1 we have 

(a) Qiita) = S'x(g«_j(to)) 

(b) q,{to) = Sy{qN_,^{ta)) 

(c) Qiita) = 5'o(g«+,(^o)) 

(d) Uto) - ^o(%+.(io)) 

(3) At time ti = T/2 for i = 0,1, ... , N/2 - 1 we have 

(a) q,{ti) = S'x((7«_,_i(ii)) 

(b) g«(ii) = ^,('Zf-.-i(ii)) 

(c) g,(<i) = 5o(9«+.(ii)) 

(d) q,{ti) = S'o(g«+,(ii)) 

(4) — {qo{t), qi{t), . . . , qpf_i{t)) is a solution of the Newton N-body problem 
forte[0,f/2]. 

Proof: If we have a trajectory with properties PI, P2, P3 then it is easy to show 
that functions qi{t) defined in P3 satisfy conditions (l)-(3). On the other hand if 
we have a functions qi{t) which have properties (l)-(3) then we define q{t) by 



From Lemma im it follows that the proof of an existence of a doubly symmetric 
choreography is equivalent to some boundary value problem for the N-body equa- 
tion. We will now formulate this problem as a zero finding problem for a suitable 
map. 

Original phase space for the planar N-body problem has AN dimensions. It 
turns out that, if initial conditions satisfy all conditions in point 2 of Lemma ITTI it 
is enough to know values of only N variables to recover rest of them. We still may 
obtain solutions of any period, to determine this we can fix the size of curve by 
fixing one variable. Hence our reduced space be {N — l)-dimensional. In the next 
paragraph we will be more specific. 

We define map E : M^-^ — > M^^, which expands, using symmetries from (2), 
initial conditions from reduced space to the full phase space: 
{xo,yo,xo,yo,xi,yi,±i,yi, . . . ,XN-i,yN-i,iN-i,yN-i)- We consider two cases: 
N ^ 4k and N ^ 4:k + 2. 
For iV = 4fc we set 



(4.2) 




□ 



fe-1 fe-1 



E{xo X x {xk,yk)) = (0, a,io,0) x J]^ (x^, y^, i^, y^) 




x{xk,0,0,yk) X Y\_i^k~i,-yk-t,-Xk~i,yk-i) 



fe-1 



x(0, -a, -io,0) X -j/i, -ii, -y^) 



i=l 

fe-1 



x{-Xk,0,0,-yk) X Y\_i-Xk-i,yk-t,ik-t,-yk~i) 
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For = 4fc + 2 we set 
fc 



E{±o X Yi^^^'yi'^^'Vi)) = iO,a,io,0) x Y\.{xi,yi, ±i,yi) 

L i=l 

fe-1 

Y\_{xk-i, -yk-i, -ik-i,yk-i) X (0, -a, -io,0) 

1=0 

fc-i 

Y\_{-Xi,-yi,-Xi,-yi) X ]^(-.Tfc_i,?;fc_i,ifc_i, 



i=l 

fe-1 

X 

1=0 

k k-1 

X 

i=l 



In both cases a is a parameter fixing a size of an orbit. 

Again, like for the Eight, for every such initial condition exists a solution of 
the A^-body problem defined on some interval. To every initial configuration, 
following that solution, we associate, if it exists, a configuration in which for the 
first time 

• for N = Ak: bodies k and fc — 1 have equal x coordinate [xk ~ Xk-i), 

• for iV = 4fc + 2: fc-th body is on the X axis {yk — 0). 

This defines the Poincare map P : M"*" D — > M^^. 

Now, we define the reduction map R : R"'^ — > M^^^ in such way that R has 
zeroes in points satisfying conditions (3a) and (3b) in Lemma l4.1l and only in such 
points. Observe that we don't have to worry about conditions (3c) and (3d) in 
Lemma 14.11 because from the properties of if follows that if (2c) and (2d) 

holds, then (3c) and (3d) are satisfied for any ti. 
For = 4fc we set 

/N-l \ 

R I n {xt,yi,Xi,yi) I = {yk +yk~i,Xk +Xk-i,yk - Vk-i) 

fe-2 

Yl^ixi - X2k-i-i,yi + y2k-i-i,Xi + X2k-i-i,yi - mk-i-i) 



fe-2 

X 
i=0 



For = 4fc + 2 we set 

(N-l \ fe-1 

{xi,yi,Xi,yi) = {vk} X Y\_i^i - X2k-i,yi + y2k~i,Xi + X2k-i,yi - mk-i) 
i=0 / 1=0 

We define map $ : R^^^ D E-^{VL) — > R^-^ by 

<^ = Ro P o E. 

Theorem 4.2. If for some x G R^^^ $(a;) = 0, then there exists a trajectory with 
properties PI, P2, PS. 

Proof: Let x be point in which we have $(a;) = 0. Then E{x) satisfies condition 
(1),(2) in Lemma [4. II and from construction of $ we obtain that there exist tra- 
jectories {qi{t)) starting with initial condition E(x) and ending with x such that 
R{x) = 0, hence condition (3) in Lemma [4. II is satisfied. Now the assertion follows 
from Lemma [4. II □ 
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Figure 6. Gerver orbit - initial position 




3 



Figure 7. Gerver orbit - final position 

5. Existence of the SuperEight - Gerver orbit 

The Gerver orbit (Figure IHIEJ is a choreography with 4 bodies forming a linear 
chain. It's the simplest trajectory after the Eight. 
With a computer assistance we proved the following 

Theorem 5.1. The Gerver SuperEight exists and is locally unique (up to obvious 
symmetries and rescaling). 

We show the existence of SuperEight using an approach described in Sectional 
i.e. we verify assumptions of Thcorcm l4.2l Below we give some details. 
We set 

(5.1) E{xi,xo, ill) = (0, a, xq, 0, xi, 0, 0, yi, 0, -a, -xo, 0, -xi, 0, 0, -yi), 

where a is a parameter fixing the size of the orbit. 
(5.2) 

R{xa,yo,xo,yo,xi,yi,xi,yi,X2,y2,i2,y2,X3,y3,X3,y3.) = {yi+yo,xi+xo,yi-yo), 
Poincare section S is defined by 

(5.3) S{xo,yo,xo,yo,xi,yi,xi,yi,X2,y2,X2,y2,X3,y3,xs,y3) = xi - xq = 0. 

At first we have proved an existence of a zero of $(a;) using an interval Newton 
method, but in the paper we present data from proof based on Krawczyk method. In 
case of Gerver orbit the choice of the method isn't so important, because if we take 
a time step small enough or smaller set [X], then the computed matrix [f§([-'^])] 
becomes invertible and the proof usually goes through. But in general the main 
problem of application of an interval Newton method is that of an invertibility of 
[^([X])] . We avoid this using the Krawczyk method. 
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Initial values 


X 


(1.382857, 1.87193510824,0.584872579881) 


a 


0.157029944461 


[X] 


x+ [-10-^10-^)3 



Table 5. Data from the proof of existence of Gerver SuperEight. 
Initial values. 



Computed values 



-2.15400 0.257911 0.786925 
-0.08163 0.293713 0.043565 
0.939059 -0.10027 0.158399 

■[-2.87020e - 09, -2.26613e - 09] 
[-1.21155e - 08, -1.06812e - 08] 
[-5.45542e - 08, -5.10016e - 08] 

6.040646 - 10" 
1.43432e-09 
3.55268e-09 



C 
diam 

lf(m) 

X(x,[X],$) 
diam K{x, [X],$) 



[-0.1664,-0.1657] 
[-0.1764,-0.1750] 
[0.87189,0.87534] 



[0.39070,0.39119] 
[3.52548,3.52654] 
[-0.0867, -0.0842] 



[0.71771,0.71790]' 
[-0.0968,-0.0964] 
[1.99599,1.99697] 



■ [1.382857036247056692, 1.382857041633411832] ' 

[1.871935113301492981,1.871935114053588922] 
[0.5848725887384301769,0.5848725902808686872] 

5.3863551396915454466 - 09' 
7.520959410811656198e- 10 
1.54243851024915557e - 09 



Table 6. Data from the proof of existence of Gerver SuperEight. 
Matrix C and results of computation. 



To check the assertion (2) in Theorem l2.2l we used the C-'^-Lohner algorithm with 
order r = 6. The time step was set to h = 0.002. As matrix C we used an inverse 

of the monodromy matrix computed nonrigorously in a point a;, i.e. C = 
Tables |S1 and El contain numerical data from this proof. 

6. Existence of the 'Linear chain' orbit for the 6 bodies 
Figure El displays a linear chain choreography with 6 bodies. 
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1 ^ \ 










^ 4 
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Figure 8. 'Linear chain' orbit for the 6 bodies - initial position 




Figure 9. 'Linear chain' orbit for the 6 bodies - final position 



In this section we report about the computer assisted proof of the following 

Theorem 6.1. The linear chain for 6 bodies exists and is locally unique (up to 
obvious symmetries and rescaling) . 

Again, we show an existence of this orbit using an approach described in Section^ 
with some minor changes. To speed up calculation and to increase accuracy we take 
into account in computation that all time q^iit) — —qo{t), qi(t) = —qi(t) and q5{t) = 
—q2{t) (we change equation 1)1. l|l doing a suitable substitution). So, our full space 
for ODE has now 12 dimensions. We use also a different time parameterization (we 
do time shift of j of the period). To use exactly approach described in Section0]we 
interchange axes (see Fig. |Sland|5J). From Lemma llTl we will get, in this coordinate 
frame, a doubly symmetric periodic solution q{t). It is easy to see that q = q{t— ^) 
is then a solution sharing needed symmetries in the original coordinate frame. All 
data are given in frame with interchanged axes. 

We set 

(6.1) E{iQ,xi,yi,xi,yi) = {Q,a,xo,0,xi,yi,xi,yi,xi, -yi, -xi,yi), 

where a is a parameter fixing the size of the orbit. 
(6.2) 

R{xo,yo,xo,yo,xi,yi,xi,yi,X2,y2,X2,y2) = (ii,xo - X2,yo + y2, io + ±2, 2/o - 2)2). 
Poincare section S is defined by 

(6.3) S{xQ,yo,xo,yo,xi,yi,xi,yi,x2,y2,i2,y2) = yi = 0- 

To find a zero of $(x) we use Krawczyk method. We check the assertion (2) 
in Theorem 12.21 using C^-Lohner algorithm (ZLoj of order r = 9 and a time step 
h = 0.0025 for a computation in point x and h = 0.001 for a computation on set [X] . 
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Initial value 



[X] 



diam [X] 



-0.635277524319 
0.140342838651 
0.797833002006 
0.100637737317 
-2.03152227864 

1.887041548253914 



-0.635277525319,-0.635277523319]' 
[0.140342837651, 0.140342839651] 
[0.797833001006, 0.797833003006] 
[0.100637736317, 0.100637738317] 
-2.03152227964, -2.03152227764] 



2.0e 


-09 


2.0e 


-09 


2.0e 


-09 


2.0e 


-09 


2.0e 


-09 



Table 7. Data from the proof of existence of linear chain for 6 
bodies. Initial values for Krawczyk method. 



As matrix C we used an inverse of the monodromy matrix computed nonrigorously 

in a point x, i.e. C = 

In Tables and ISI we present data from computation of Krawczyk method. 

7. Some technical data 

All computations were performed on AMD Athlon 1700XP with 256 MB DDRAM 
memory, with Windows 98SE operating system. We used CAPD package .Capdi and 
Borland C++ 5.02 compiler. 

In the listing below r is an order and h is a time step used in the C^-Lohner 
algorithm |ZLo| . 

The computation times for The Eight, h — 0.01, r = 7 

• in point x : 1.417 sec 

• for set [X] : 2.66 sec 

• convexity : 1.15501 sec 

For the proof of an existence of Gerver solution in 4-body problem we used r — 6, 
h — 0.002. The computation times for both x and set [X] were approximately equal 
to 30.5 seconds. 

For the proof of linear chain for 6-body problem 

• computation of Poincare Map for set [X] took 57.5 seconds with h = 0.001 
and r = 9 

• computation for x took 23.8 seconds with h — 0.0025 and r — 9 
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Computed value 



diam $(a;) 



diam K{x,[X],^) 



[-3.13119579096586- 
-4.528821762050939e 
-1.063704679893362e 
-3.0841051934515926 
-1.2037260077590876 



11,3.1562770627005856- 11] ' 
~ 12,4.5747572396948046- 12] 

- 11, 1.0514700221619936- 11] 

- 11,3.1174951509171936- 11] 

- 11, 1.193112275643671e- 11] 



6.2874728536663866- 11 
9.1035790017457436- 12 
2.1151747020553566- 11 
6.2016003443687856- 11 
2.3968382834027586- 11 



[-0.6352775243616679557,-0.6352775242763283314]' 
[0.1403428386430521646,0.1403428386590999943] 
[0.797833001999263769,0.797833002012834469] 
[0.10063773728817425324,0.1006377373457752189] 
[-2.031522278710178764,-2.031522278575771612] 



8.533962425616437036- 11 
1.6047829731746787736- 11 
1.3570700119203138466- 11 
5.7600965663873182626- 11 
1.3440715207480025136- 10 



Table 8. Data from the proof of an existence of linear chain for 
6 bodies. Results of computation of Krawczyk method. 



The programm performing the proofs is available on http://w'ww. ap.krakow.pl/ ~tkapela 
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